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Abstract 

Galaxy clusters arc promising targets for indirect dark matter searches. Gamma-ray signatures 
from the decay or annihilation of dark matter particles inside these clusters could be observable 
with the Fermi Large Area Telescope (LAT). Based on three years of Fermi LAT gamma-ray 
data, we analyze the flux coming from eight nearby clusters individually as well as in a 
combined likelihood analysis. Concentrating mostly on signals from dark matter decay, we 
take into account uncertainties of the cluster masses as determined by X-ray observations 
and model the cluster emission as extended sources. Searching for different hadronic and 
leptonic decay and annihilation spectra, we do not find significant emission from any of the 
considered clusters and present limits on the dark matter lifetime and annihilation cross- 
section. We compare our lifetime limits derived from cluster observations with the limits that 
can be obtained from the extragalactic gamma-ray background (EGBG), and find that in 
case of hadronic decay the cluster limits become competitive at dark matter masses below 
a few hundred GeV. In case of leptonic decay, however, galaxy cluster limits are stronger 
than the limits from the EGBG over the full considered mass range. Finally, we show that 
in presence of dark matter substructures down to 10~^ solar masses the limits on the dark 
matter annihilation cross-section could improve by a factor of a few hundred, possibly going 
down to the thermal cross-section of 3 x 10~^^ cm'^s"^ for dark matter masses < 150 GeV 
and annihilation into bb. As a direct application of our results, we derive limits on the lifetime 
of gravitino dark matter in scenarios with i?-parity violation. Implications of these limits for 
the possible observation of long-lived superparticles at the LHC are discussed. 
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weniger@mppmu.mpg.de (C. Weniger) 



1 Introduction 



Galaxy clusters are the most massive gravitationally collapsed objects in the Universe. Histor- 
ically, the kinematical study of the Coma cluster provided the first indication for the existence 
of dark matter Q. Further gravitational evidenc; for dark matter was found at many places, 
e.g. in stellar rotation curves of spiral galaxies or as baryon acoustic oscillations imprinted in the 
cosmic microwave background (for reviews on particle dark matter see Refs. 0-0] )• By now, the 
ACDM scenario is the standard framework for cosmology, leaving open the question of what the 
nature of the dark matter particles is. Theoretical models for dark matter predict a large variety 
of possible non-gravitational signatures that would help to reveal its properties. However, despite 
lots of efforts, none of these signals has been unambiguously detected so far. 

If dark matter is made out of WIMPs (Weakly Interacting Massive Particles), their efficient 
self-annihilation in the early Universe would explain the observed dark matter density. Today, the 
same annihilation process could contribute to the measured cosmic-ray fluxes; a clear detection 
of the annihilation products would reveal information about the dark matter particle's mass and 
interactions. Similar signatures could be produced if dark matter is unstable and decays, providing 
us with information on the lifetime of the dark matter particle. Indirect searches for dark matter 
are aiming at seeing such annihilation or decay signals above the astrophysical backgrounds. 
These searches typically concentrate on photons or neutrinos, which carry spatial and spectral 
information about their origin, and on different anti-matter species with their relatively low 
astrophysical back grou nds. Together with satellite galaxies of the Milky Way js-ll] and the 
Galactic center l^- 3| , nearby massive galaxy clusters are among the most promising targets for 
indirect dark matter searches by means of gamma rays 

A theoretically well motivated example for decaying dark matter is the gravitino '4'3/2j which 
appears in locally super symmetric extensions of the Standard Model. In scenarios where -R-parity 
is mildly violated and the gravitino is the lightest superparticle (LSP), thermal leptogenesis. 



gravitino dark matter and primordial nucleosynthesis are naturally consistent [3ll |. Within this 
framework, the gravitino would decay with cosm ologica l lifetimes [32] , making its decay products 
potentially observable in the cosmic-ray fluxes j33l444i| . For gravitino masses < 100 GeV, the 
most prominent feature in the decay spectrum is an inten se g amma-ray line, produced by the 



two-body decay into neutrinos and photons, '03/2 ~^ 7^ [SJ]- Dedicated searches for such a 
feature in the current gamma-ray observations of the Fermi LAT [45[ exist in the literature, see 
Refs. MA 

and their null results were used to put lower limits on the gravitino lifetime around 
6 X 10^^ s [471]. However, for larger gravitino masses > 100 GeV, the branching ratio into gamma- 
ray lines is strongly suppressed, and instead decay modes like "03/2 ~^ W^t^ and 1/^3/2 — Z^v 
produce a gamma-ray flux with a broad continuous energy spectrum. It is this flux that could 
potentially show up in observations of galaxy clusters, whereas the observation of the gamma-ray 
line in galaxy clusters would be difficult due to the limited statistics. 

In general, dark matter lifetimes of the order of 10^^-lO^^s, which is in the ballpark of what 
is accessible experimentally, are obtained when the symmetry responsible for the dark matter 
stability is violated by dimension six operators generated close to the grand unification scale [48- 
501 ]. Indeed, models of this kind were proposed to explain the "excesses" observed by the 



PAMELA [5l|], Fermi LAT 



531 ] and H.E.S.S. [5^, |55l] experiments. To avoid the stringent 
anti-proton limits 5l|, l56|, the decay should be mostly leptophilic (see e.g. Refs. ST-G^]), and 
one typical decay mode that could well reproduce the locally observed fluxes is the decay into 
muons, ijj — )• n~ , with a large dark matter mass around ~ 3 TeV and lifetimes around 
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2 X 10^^ s 6^, |65[. Inside galaxy clusters, due to inverse Compton scattering (ICS) on the cosmic 
microwave background (CMB), almost all of the kinetic energy of the produced high-energetic 
electrons and positrons is transferred into gamma rays, with energies up to 0(100 GeV). This 
makes possible the investigation of the decaying and annihilating dark matter interpretations of 
the excesses by galaxy cluster observations. 

A dedicated search for dark matter annihilation signals from galaxy clusters, using 11 months 
of Fermi L AT data, was carried out in Ref . [2^ . The null result of this search was used to derive 



limits on the dark matter annihilation rate into bb and ii~ ^'^ . In Ref. [271] these results were 
translated into limits on the dark matter decay rate, and it was demonstrated that galaxy cluster 
observations give strong constraints on the dark matter lifetime, superior to the limits that could 
be obtained from satellite galaxy observations, and of the order of the limits that can be derived 
from the extragalactic gamma-ray background. This makes galaxy clusters promising targets 
when searching for signals from dark matter decay. Concerning WIMP dark matter, taking into 
account the expected boost of the annihilation signal due to dark matter substructures in the 
cluster halo, limits can potentially go down to the cross-section expected from thermal freeze- 

uster emission as seen by the Fermi LAT, H.E.S.S. 
281, some applications to annihilating and decaying 
66 -B. 



out [29|, |30(] . Further studies of the galaxy c 



and MAGIC were presented in Ref. [23|, [24, 



dark matter models were discussed in Refs. 

Besides the large amount of dark matter, it is known from radio observations that galaxy 
clusters are also a host for energetic cosmic rays, which can be accelerated during the process of 
cluster formation by mergers or accretion shocks. Proton-proton collision as well as the ICS of an 



energetic electron population can produce a possibly observable gamma-ray flux (see e.g. Refs. j26l . 



30[). Such a flux should however be finally distinguishable from a dark matter signal through the 



analysis of the energy spectra if the statistics is high enough 



In this paper, we analyze the gamma-ray flux from eight galaxy clusters as measured by 
the Fermi LAT since Aug 2008, and we present constraints on the dark matter lifetime and 
annihilation cross-section. We analyze the different target clusters individually as well as in a 
combined likelihood approach, and search for significant gamma-ray emission as an indication for 
decaying or annihilating dark matter. In Ref. 2^, 27 1 the dark matter signal was approximated 



to be point-source like. Importantly, this approximation becomes problematic in the case of 
dark matter decay or substructure-boosted annihilation, since the extend of the expected signal 
starts to exceed the angular resolution of the Fermi LAT. To account for this, we model the 
dark matter emission as an extended source. In absence of a clear signal, we derive limits on 
the dark matter lifetime and annihilation cross-section as function of the dark matter mass, for 
different hadronic and leptonic final states. Cluster masses and the expected decay or annihilation 



signals are derived from the extended HIFLUGCS catalog [69|, l70|] which is based on ROSAT 



PSPC X-ray observations [7l|], and the corresponding uncertainties are consistently taken into 
account. We compare the obtained lifetime limits with the limits that can be derived from 
the extragalactic gamma-ray background, and we discuss the implications of our limits on the 
decaying or annihilating dark matter interpretation of the excesses. Furthermore, we will 
illustrate how the limits improve when a boost of the annihilating signal due to substructures is 
included. Finally, we apply our findings to the scenario of decaying gravitino dark matter, and 
derive new constraints on the gravitino lifetime for masses above about 100 GeV. We comment 
on implications for the possible observation of long-lived superparticles at the LHC. 
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Cluster R.A. Dec. z J^''^- J^^- 

[lOiSGeVcm-2] [lO^^GeV^cm^^i 



Fornax 


54.67 


-35.31 


0.0046 




8 Q+2.0 
0-0-2.8 


0.44 


Coma 


194.95 


27.94 


0.0232 


10.7lij 


1 q + 0.20 

l--"-'-0.31 


0.23 


A1367 
A1060 


176.19 
159.18 


19.70 
-27.52 


0.0216 
0.0114 


10.2^1° 


1 4+0.15 

9 9+0.38 
^•^-0.69 


0.23 
0.24 


AWM7 
S636 

NGC4636 
NGC5813 


43.62 
157.52 
190.71 
225.30 


41.58 
-35.31 
2.69 
1.70 


0.0172 
0.0116 
0.0037 
0.0064 


q q+1-9 
6.8l|-^ 

c 1 +0.80 

6.0lt^ 


■•-•"-0.56 
1 tc+0.29 
■•-•^-0.34 
q r+0.39 
•^•^-0.85 
9 9+1.4 
^•^-1.4 


0.22 
0.18 
0.19 
0.18 



+0.07 
-0.11 
+0.02 
-0.04 
+0.02 
-0.04 
+0.03 
-0.06 
+0.03 
-0.06 
+0.03 
-0.03 
+0.02 
-0.04 
+0.08 
-0.10 



Table 1: Galaxy clusters considered in this work, with their coordinates (equatorial J 2000.0) and 
redshift z from Ref. [fe^J . We show the integrated surface densities J^q = J^^ J''™' (f2) and J^q' = 
Iaq J™"' of the dark matter signal, obtained inside a region of 1° radius around the cluster center, 
as well as the projected angle 6s of the scale radius Ts of the adopted NFW profile. Central values 
and errors for these parameters are derived from the cluster masses in Ref [/dj . 

The rest of the paper is organized as fohows: In section 2 we present our galaxy cluster analysis 
of the Fermi LAT data. We discuss the expected dark matter signals, our treatment of the LAT 
data and the details of our statistical analysis. In section 3, we shortly review how limits on the 
dark matter lifetime from the extragalactic gamma-ray background are obtained. Our results 
and their discussion are presented in section 4. Finally, section 5 is devoted to gravitino dark 
matter, where we present limits on the gravitino lifetime for gravitino masses above 100 GeV, 
as well as the implied limits on the decay lengths of next-to-lightest superparticles (NLSP) at 
particle colliders. We conclude in section 6. 

Throughout this work we assume a ACDM cosmology with the parameters f^A = 0.728 and 



h = i/o/lOOkm s-i Mpc"^ = 0.704 [T^. 



2 Galaxy Cluster Analysis 

The eight galaxy clusters that we consider in this work are summarized in Tab. [TJ They are 



selected from the extended HIFLUGCS X-ray catalog [69|, |70|] in order to yield large signals 
from dark matter decay, but are also among the best clusters when searching for signals from 
dark matter annihilation. Galaxy clusters with potentially large signals that we disregard are: 
Ophiuchus, A3627 and 3C129 because they lie too close to the Galactic plane, Centaurus, M49 
and A2877 because of issues with our adopted model for Galactic diffuse emission^ and Virgo 



(M87) and Perseus because of the presence of bright gamma-ray sources at their center |73l. l74|. 



^The positions of the Centaurus and M49 clusters unfortunately coincide with sharp edges in our Galactic 
diffuse emission model, gal_2yearp7v6_v0; the region near A2877 contains a large number of faint sources that are 
not part of the 2FGL. In all three cases the background fits are unreliable, and we neglect these targets from our 
analysis. Including them would improve our overall limits. 
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2.1 Dark Matter Signal 

The gamma-ray flux from dark matter annihilation or decay that is expected to be seen in galaxy 
cluster observations factorizes into an astrophysical part, which contains information about the 
dark matter distribution pdmi and a particle-physics part, which is universal for all observed tar- 
gets. Assuming a spherical dark matter halo, the astrophysical factor, jdcc./ann. ^^-j^ j^g^. (jgpgj^ds 
on the cluster-centric angle and is given by a line-of-sight integral. In the case of dark matter 
decay, the signal flux reads 

''^^^W = 7r::^^ / dsp,^is,n), (i) 

/l.o.s. 



dE dO, 47rm^ dE 



while in the annihilation case, it is given by 



dJsig , , (av) dN^/ 



dE dn^ ' Svr m\ dE „ 



W = ^^-7T?/ dspi^{s,^). (2) 



(6») 

Here, denotes the dark matter mass, while and {av) are the dark matter lifetime and total 
annihilation cross section, respectively. The energy spectrum of gamma rays produced in the 
decay/annihilation is given by dN^/dE. Note that the energy spectrum dN^/dE = dN^^^^/dE-\- 
dN}^ /dE includes prompt gamma rays that are directly produced in the decay or annihilation 
process (final-state radiation, vr'^ — ?• 77 etc.) as well as the gamma rays that originate from ICS 
losses of from dark matter on the intra-cluster radiation field. We calculated the energy spectra 
of gamma rays and electrons with the event generator Pythia 6.4.19 [tH], and cross-checked our 
results with the analytic expressions presented in Ref. [76[. 

Inverse Compton scattering. Electrons and positrons produced in the decay or annihilation 
of dark matter particles inside galaxy clusters suffer inverse Compton and synchrotron losses when 
interacting with the intra-cluster radiation field. The dominant component of this radiation field 
is in most cases the CMB; other contributions, which can become relevant close to the cluster 
center, are the starlight, dust radiation and the intra-cluster magnetic field (see discussion in 
Ref. [30|). In case of dark matter decay or substructure-boosted annihilation, the possible impact 
of these additional components on our results is small, as we will exemplify below for the Coma 
cluster; if not stated otherwise, we consider the CMB only throughout this work. 

The average energy spectrum of gamma-rays with energy E^ that are generated by the inverse 
Compton scattering of one electron with an initial energy of Eq is given by 



dE^ Jq J^^ dEy bioss{Ee) 

Here, /cMB(e) is the CMB energy spectrum with temperature Tcmb = 2.725 K, and da^'^ /dE^ 
denotes the differential cross section of inverse Compton scattering of an electron with energy Ee 
when a CMB photon with energy e is up-scattered to energies between E^ and E^ + dE^. Due 
to the very low energy of the CMB photons, the center-of-mass energy of the processes we are 
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interested in is always smaller than the electron mass, which allows us to use the non-relativistic 
limit of the Klein-Nishina equation in our calculation^ (see e.g. Ref. 77]): 



where <tt = 0.67 barn is the Compton scattering cross section in the Thomson limit, je = E^/me 
is the Lorentz factor of the electron, mg = 511 keV is the electron mass, and q = E^/E^'^^ with 
^max ^ 47ge. Eq. @ holds in the range e < Ey < E™^^, in the limit where down-scattering is 
neglected. 

In the above equation, bioss{Ec) is the energy loss rate of an electron with energy E^. Typically, 
the main contribution comes from ICS on the CMB, but we can also include synchrotron losses 
on the intra-cluster magnetic field; 6ioss = &ics + ^syn- In the non-relativistic limit, the ICS and 
synchrotron losses read 

blCs{Ee) = -O-Tle / de e/cMB(e) and 6syn(£^e) = 3'^T7e^ , (5) 

^2 . ' 

=PCMB 

respectively. In order for the magnetic field B to dominate the CMB energy density pcMB (namely, 
B'^/2 > pcmb) in galaxy clusters it has to exceed the critical value -Bcmb = 3.2^G (assuming 
redshifts z <^ 1). The energy spectrum of the dark matter induced ICS radiation would then 
scale like oc {1 + (5/3.2/zB)2)-i. 

The energy loss time rioss = Ee/Eg of electrons with 100 GeV-10 TeV energies, as relevant 
for our work, lies in the range of 10~^-10~^ Gy. This is much shorter than the cosmic-ray 
relaxation times in galaxy clusters which typically are of order 1-10 Gy 3C. 7^. The propagation 



scale corresponding to > 100 GeV electrons is expected to be < 1 kpc [79|, l80|], which for Mpc 



distances is well below the angular resolution of the LAT [811]. We hence neglect effects of cosmic- 
ray transport and consider that are loosing all their energy instantaneously where they are 
produced. In this limit, the angular profile of the ICS signal is identical to the angular profile of 
the prompt radiation. This is opposite to the case of dwarf galaxies, where due to their proximity 
propagation effects have to be taken into account in general, see e.g. Ref. 0]. 

After some algebra, one can finally show that the energy spectrum of ICS radiation emitted 
from a single electron with an initial energy Ef, is given by the expression {B — t- 0) 



dN}p 9me /"i , 1 r 92 23/2 2 4 A ] /cmb (4^ 

^9 ^ <^ 77^ - 779^/^ - 7779^ + t;Q ^ - Ti' log q \ . (6) 



47e 



dE^ 327I 7 1 g5/2 \525 3 25 7 5 j pcMB 



A subsequent convolution with the energy spectrum dNf>/dEf^ of electrons and positrons yields 
then the ICS contribution to the gamma-ray spectrum, dN}p/dE, in Eqs. ([T]) and ([2|). 

2.2 Dark Matter Distribution 

As discussed above, the expected dark matter signals depend crucially on the dark matter profiles 
p(r)dm of the target galaxy clusters. We assume throughout this work that the smooth component 



This approximation breaks down for electron energies above m^/TcMB ~ 10^ TeV. 
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of the dark matter halo follows a Navarro- Frenk- White (NFW) profile 46|, [8 



Pdm{r) = ^ , (7) 

r/rs{l + r/rsY 

where the scale radius and the density normalization have to be determined from observa- 
tionsH The cluster mass Ma inside a cluster-centric radius ta is defined such that the average 
density inside ta equals A times the critical density of the Universe, pc (typically A ~ 100-500). 
To determine the parameters of the NFW profile from Ma and ta, we adopt the observationally 
obtained concentration-mass relation from Ref. [Sf 



where the concentration parameter c = r^ir/vs relates the virial radius r^ir as defined by A = 
Avir — 98 (see appendix of Ref. [s^ and references therein) to the scaling radius , and Mq — 
2.0 X 10^^ kg denotes the solar mass0 The virial radius r^ir is then related to ta via 



/(r./rA) = -^firs/r^ir) , (9) 

^vir 

where 

f{x) = [ln{l+x~^) - + , (10) 

and 

With Eqs. ([8|), Q and ([TT]) . one can find and ps as a function of M/\ and r/\. 

Using the cluster masses as derived from ROSAT PSPC X-ray observations in the extended 
HIFLUGCS catalog [t^], we calculate the signal surface densities J^""^- and J'^'^'^' as described 
abovell Our results, as function of the cluster-centric angle 9, are plotted in Fig. [TJ These 
profiles are used to model the extended dark matter signal in our analysis. For convenience 
and comparison with previous work, we show results for J^^'^- and J'^'^'^- integrated over a cluster- 
centric region of 1° radius in Tab.[TJ There, we also indicate the projected scaling angle Og = Vg/D, 
where the distance to the cluster is given by D ~ zc/Hq. The signal uncertainties shown in Tab.[l] 
are directly derived from the mass uncertainties in Ref. [t^] and as large as a factor of two in some 
cases. Within the error bars our results agree largely with what was found in Refs. 2^, 27] based 



on the initial HIFLUGCS catalog 1691]. Besides the uncertainties from the fits to the X-ray profiles 
that were already discussed in Ref. |69l]. the mass ranges quoted in Ref. [t^] take additionally into 
account uncertainties in the X-ray temperature profile, which leads to somewhat larger error-bars 
compared to Ref. [g^. 



■^We find that using an Einasto profile [83l . |8J] with similar M200 and r2oo (with a = 0.17, r_2 = Ts and 
p_2 — Ps/4.2 [s^l) leaves the resuhs for decaying dark matter essentially unchanged, whereas the fluxes from dark 
matter annihilation as summarized in Tab. [T] are increased by about ~ 30%. 

^Varying the prefactor in Cvir in the range 9±2 leaves the fluxes from decaying dark matter essentially unchanged, 
whereas the annihilation fluxes increase or decrease by about 30-40%. Note that CDM simulations favor a somewhat 
smaller concentration of about ~ 6 at Mvir = 10"/i"^Mq [ll,!!!! . 

®The values of M500 and rsoo are rescaled to our adopted Hubble constant, and M500 is reduced by the gas 
fraction /gas indicated in Ref. [t^I ■ We checked that using M200 — \/5/2 M500 and r2oo — 1/6/2 rsoo as a starting 
point would increase the predicted fluxes by 10-20%. 
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Figure 1: Left panel: profile of the dark matter decay signal as function of the cluster-centric angle 
9 (PSF effects not included). The dotted line shows the isotropic Galactic contribution to the dark 
matter signal. Right panel: the same for the dark matter annihilation signal. Solid lines show 
the signal coming from the smooth halo component alone, dashed lines include effects from dark 
matter substructures, which boosts the signal at angles around 1° (as discussed in Section \2.3( l. For 
comparison: the angular resolution of the LAT (P7S0URCEJ/6) in terms of the 68% containment angle 
is 6° at 100 MeV, 0.9° at 1 GeV and 0.2° at 100 GeV 0/. 



2.3 Signal Boost from Dark Matter Substructures 



A prediction of the cold dark matter paradigm is the hierarchical structuring of dark matter halos. 
Dark matter substructures inside of galaxy cluster halos are observationally known to exist down 
to the scale of dwarf galaxies, lO^M©; for thermally produced WIMPs they are predicted to 
continue down to free streaming masses of about W~^Mq and below [s^, [oO] (for a discussion 
of possible ranges depending on the dark matter model see Ref. [9l[). Since the dark matter 



annihilation signal depends on the dark matter density squared, the existence of substructures 
can boost the annihilation signal considerably with respect to the signal from the smooth halo; 
the details depend on the mass function of substructures, the concentration mass relation and 
the radial distribution (see e.g. Ref. [30(] for a recent discussion). Dynamical friction and tidal 
stripping near the cluster center lead to a local depletion of substructures that results in a relative 
enhancement of the boosted signal in the outskirts of the main halo. In general, the boosted 
signal is expected to be considerably more extended than the signal coming from the smooth 
main halo alone. Deriving the magnitude of the signal boost relies on extrapolations of numerical 
simulations for dissipationless DM 92H94i| over many orders of magnitude in the substructure 
mass. In the literature, predictions for the substructure boost have not yet converged; in case of 

0, Q up to ~ 1000 0, [sjl were recently 



galaxy clusters, signal boosts in the range of ~ 10-50 
discussed. As mentioned above, the actual values strongly depend on the adopted subhalo mass 
fraction (which is partially correlated with the value of erg used in the underlying simulations), 
the subhalo mass distribution functions and the adopted halo concentration. Furthermore, for 
individual clusters, precise predictions appear to be difficult since the halo-to-halo scattering of 
the substructure fraction, which roughly correlates with the concentration of the particular halo, 
can be quite large and 0(1) jsH]- 

In the present paper, we adopt an optimistic scenario and estimate the boosted dark matter 
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signal following Ref. [30|]: Based on the high-resolution dissipationless dark matter simulations 



of the Aquarius project 92] (which features a realtively large subhalo fraction), the boost of 



the dark matter annihilation signal was determined in Ref. 93[| for a Milky Way sized halo. 
Extrapolating the mass of the smallest subhalos down to Mn^ = 10~^Mq, an increase in the 
overall luminosity L = f dVp"^^ of about 230 was found, being mostly due to a signal enhancement 
at large galactocentric distances. The luminosity due to substructures inside a radius r is well 
fitted by [iO] 

^sub(< r) = aoC(M2oo) ^200sm(M2oo) , (12) 

fix) = aix'^^ , 



C(M2oo) 



0.023 — ^ 



where x = r/r2oo, oq = 0.76, ai = 0.95, 02 = —0.27 and ac = 0.226. Here, L200sm denotes the 
luminosity of the smooth halo component inside r2oo alone. The only free parameter is the cutoff 
scale for the dark matter subhalo mass, which we fix to = 1O~^M0. The parameters M200 
and r2oo are directly determined from the adopted NFW profile for each cluster. The overall 
boosted dark matter signal can then be calculated from Eq. (jl2p . 

Our results for the signal profile from dark matter annihilation in presence of substructures are 
shown in the right panel of Fig. [T] by dashed lines, the signal from the smooth halo alone is shown 
by solid lines. In presence of substructures, the annihilation signal extends to radii of around 
1°, below 0.01° it is still dominated by the smooth dark matter halo. The boost factors that we 
obtain for the different considered galaxy clusters inside an opening angle of 6*200 = f2m/D are 



in the range 500-1200, in agreement with Ref. |30l |. Similar large values were recently also found 



in Ref. [96( (however, see discussion above). For different values of the cutoff M\\^ the boosted 

-o.r" 

lim 



signal scales like oc M,- 'J'^^^, whereas its angular profile remains unchanged. 



2.4 Data Analysis 

The gamma-ray events entering our analysis are selected from the P7S0URCE_V6 event class of 
the Fermi LAT data measured between 4 Aug 2008 and 21 Jul 2OI1H From all events recorded 
by the Fermi LAT, we select those with energies between 400 MeV and 100 GeV and apply the 
zenith angle criterion 9 < 100° in order to avoid contamination by the Earth's AlbedoQ For each 
galaxy cluster, we consider photons events in a 10° x 10° squared region centered on the cluster 
position. These events are binned into a cube of 0.1° x 0.1° pixels with 24 logarithmic energy 
bins. The lower end of the considered energy range is somewhat larger than what was used in 
previous works, e.g. Refs. j^, [2^: Below energies of 400 MeV, the point spread function (PSF) of 



the LAT becomes of the size of our considered target regions 8l[, whereas in the considered dark 
matter scenarios no relevant gamma-ray fiuxes below ~ 400 MeV are expected; this motivates 
our choice. 

For the diffuse background fluxes we take the isotropic emission and the galactic foreground 
model templates currently advocated by the Fermi LAT collaboration for point source analysis 



^The event data as well as the corresponding information about the instrument response functions P7V6 can be 
obtained from http: //f ermi . gsf c .nasa. gov/ssc/data/ We checked that using the event class P7CLEAN_6 instead 
leads to results that are similar to what is presented in this paper. 

^These selections are made using the Fermi Science Tools v9r23pl. For the cuts in gtmktime we took DATA_QUAL==1 
as well as the RIO-based zenith angle cut. 
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(iso_p7v6source and gal_2yearp7v6_v0). The galactic foreground model contains several ded- 
icated spatial templates to model diffuse emission that is not accounted for by the GALPROP 
code {e.g. for Loop I and the Galactic Lobes [o^DII These spatial templates exhibit sharp edges, 
and we exclude clusters that coincide with these edges (Centaurus and M49) from our analysis 
in order to avoid a bias of our results. On top of the diffuse templates, we add the point sources 
from the second Fermi LAT catalog 2FGL [98(. We include all point sources within a radius of 
12° around the cluster centers. Some of the sources lie outside of our target regions but might still 
contribute due to the large point-spread function (PSF) of the Fermi LAT at low energies. We 
furthermore (and conservatively for the purpose of deriving limits) assume that gamma-ray emis- 
sion due to e.g. shock accelerated cosmic rays inside the cluster is absent and attribute possible 
observed fluxes entirely to dark matter. 

We use the profile likelihood method to fit the data and derive limits (o^ . 10Cl( |. The corre- 



sponding likelihood function C is — for an individual clusters j — given by Cj{c^fj,^) = niP(c^|//j), 
where -P(c|/x) denotes the Poisson probability to observe c events when n are expected. The num- 
ber of expected counts /i^ = /^^(q:) that is predicted for an energy /spatial bin i is a function 
of the model parameters ex. These numbers are in principle obtained by a convolution of the 
above model fluxes with the instrument response function of the Fermi LAT. In this work, the 
convolution with the PSF is done using gtsrcmaps from the Fermi Science Tools. Like in most 
of the existing analyses of Fermi LAT data (for exceptions see Refs. 0, 0, |43|) we will neglect 
the small but finite energy dispersion of the LAT, which would not significantly affect the broad 
energy spectra that we are considering. Finally, the best-fit model parameters for cluster j are 
obtained by maximizing Cj{c^\fj,^ {a)) with respect to the model parameters a. 

As a first step, we fit the data extracted from our eight target regions with the background 
model only. The free parameters in the fit are the normalizations of the two diffuse background 
templates, as well as the normalization and spectral index of all point sources inside a 5° radius 
around the cluster position or with TS- values larger than 9. In Fig. [2] we show for two exemplary 
clusters the residual maps that we obtain after subtracting our best-fit background models from 
the data, integrated over all energies; the corresponding energy spectra of the individual back- 
ground components are shown in Fig. [3j For other clusters, we obtain similar results. The figures 
indicate that the adopted background models are sufficient to model the observations. 

We then include the potential dark matter signals in the fits. The individual cluster signals 
are modeled as extended sources; their surface densities follow from Eqs. ([1]) and ([2]) and are 
plotted in Fig. [T] We neglect the smooth Galactic and extragalactic components {cp. Fig. [H 
dotted line) of the dark matter signal and assume that they are already accounted for by our 
two diffuse templates!! Uncertainties of the cluster masses as given in Ref. 70] translate into 



uncertainties on the integrated signal from each cluster, and into uncertainties on its angular 
shape as parametrized by 9s', the resulting errors are indicated in Tab. [TJ We find that variations 
in 9s have much less impact on our limits than variations of the integrated signal (less than 10% 
for the ranges given in Tab.[l]); for simplicity we will keep 9s at its central value when performing 
fits to the data. 

We include uncertainties of the cluster mass as a systematic error into the profile likelihood 



*See http: //fermi .gsfc .nasa.gov/ssc/data/access/lat/Model_details/Pass7_galacti c .html I 

^This is a realistic assumption, since these contributions to the dark matter signal could be ea sily m istaken as 
part of the extragalactic gamma-ray background as determined by the Fermi LAT collaboration [lOll ]. Since we 
already include a template for this extragalactic flux, a further inclusion would lead to double counting. 
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Figure 2: Residual maps after subtraction of our best-fit background models, in units of la standard 
deviations, for the case of the Fornax and Coma cluster. The maps span a 10° x 10° region and are 
centered on the cluster position, pixels are resampled to 0.5° x 0.5°, counts are summed over the full 
energy range 400 MeV-100 GeV. The count number per resampled pixel ranges between 11 and 273 
with an average of 28.0. 



method by substituting the hkelihood function of an individual cluster Cj (with a dark matter 
signal modeled according to the central values of Tab. [1]) with a likelihood function that takes 
into account the corresponding uncertainties of the integrated dark matter signal [ll| 



CjiaBu) ^ 4("dm) = max , Cj ( odm^ 1 J^f^\Jin) ■ (13) 



•Jaq \ pi-' 



An, 

Here, odm denotes the normalization of the dark matter signal (being related to the dark matter 
lifetime or annihilation cross-section), is the integrated surface density, J^q its central value 
as given in Tab. [H and ^f'^^ is the likelihood function of J^^^ for cluster j. In this work, 
we approximate Cj^^^ by a log-normal distribution that is defined according to the error bars 

in Tab. [Din order to model the uncertainties of J^^^. To this end, we fix £,^^ such that its 
cumulative distribution function equals 0.16 and 0.84 at the lower and upper errors given in 
Tab. m respectively!^ 

In our signal-l-background fit, we fix most of the 2FGL sources to their values from the 
above background-only fit; exceptions are the sources 2FGL J1037. 5-2820, 2FGL J0334. 3-3728 
and 2FGL J1505. 1+0324, which lie close to the A1060, Fornax and NGC5813 cluster positions, 
respectively. This leaves us for most of the clusters with three free parameters: the signal nor- 
malization odm and the two normalizations of the diffuse backgrounds. We checked that leaving 
more 2FGL source parameters free in the fits does not change our results significantly, but in- 
creases the computational time considerably. We scan the likelihood function ^f^^ as function 
of aDM while refitting the remaining free parameters. Upper limits at the 95% C.L. (99.7% C.L.) 
on the dark matter signal can be derived by increasing the signal until — 21og£j^*^ increases by 
2.71 (7.55) from its best-fit value. The significance of a signal can be obtained by comparing the 
likelihood values that are obtained with and without a dark matter signal. 



^"Redoing the error analysis that was performed in Ref. [7y], we found that a log- normal function describes well 
the posterior probability distribution function (pdf) of the cluster masses (as long as the polytropic index is 7 ~ 1). 
Since J^Q oc (M500)" with a ~ 1, the same holds for the posterior pdf of J^q. Assuming a flat prior for the cluster 
masses, this motivates us to adopt a log-normal function also for the likelihood function of J^q. 
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Figure 3: Energy spectra of different background components compared with data in the energy- 
range 100 MeV-100 GeV, for Fornax (left panel) and Coma cluster (right panel). Blue and red lines 
correspond to our two diffuse templates gal_2yearp7v6_v02 and isojp7v6source, respectively, the 
lower black lines show the contribution from different point sources of the 2FGL. When performing 
the fits, we only use data down to 400 MeV, but the background models continuously connect also 
to data at lower energies, as shown in the plot. 

Finally, to combine the statistical power of the different target regions and to reduce the 
impact of the cluster mass uncertainties, we performed a combined likelihood analysis of all eight 
clusters simultaneously. In this case, the combined likelihood function Ccomb is defined as the 
product of the individual likelihood functions, >Ccomb = ^j^ji where j runs over the different 
galaxy clusters. The only parameter that is bound to be identical for all targets is the dark 
matter lifetime or annihilation cross-section!^ 

Note that we use our own software to profile over the combined likelihood function in pres- 
ence of cluster mass uncertainties. These scanning routines were implemented on top of the 
Fermi Science Tools (and are independent of the routines used in the combined dwarf analysis of 
Ref. HH). 

3 Limits from the EGBG 

In case of dark matter decay, an important contribution to the gamma-ray signal always comes 
from our own Galaxy. Assuming an NFW profile (r<j = 20 kpc and pq = 0.4 GeV/ cm'^), we 
obtain J|fe|>ioo = 2.1 • 10^2 GeV cm ^ sr when averaging over the whole sky excluding the 
Galactic disk, and J^^'^gQo = 1.1 • 10^^ GeV cm"^ sr~^ at the Galactic anti-center (the maximal 

^^Note that the angular distance between the targets A1060 and S636 is only 8° and hence their target regions, 
but not the signal regions, overlap to a certain degree. We checked that when profiling over the signal normalization 
up to the 2(7 limits (in order to obtain Cj) the background normalizations are only affected at the < 1% level, 
hence the limits on A1060 and S636 remain practically statistically decoupled and the combined likelihood analysis 
is applicable. 
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isotropic component of the Galactic flux). In Fig. [T] we compare the angular profiles of the cluster 
decay signal with the contribution from our Galaxy in anti-center direction. As evident from this 
plot, the Galactic component dominates the signal already at a distance above ~ 0.5° from the 
cluster center. 

For comparison with our galaxy cluster limits, we will derive additional limits on decaying dark 
matter by requiring that the isotropic component of the Galactic signal plus the spatially averaged 
extragalactic signal does not overshoot the extragalactic gam ma-ray b ackground (EGBG) as 
derived by the Fermi LAT collaboration 10l| (see e.g. Refs. 33, 1021 - lod ]). In the calculation of 
the prompt signal component, we fully take into account the Galactic (in anti-center direction) 
and the red-shifted extragalactic sign al flux, and we employ for completeness the inter-galactic 
background light model of Ref. 107l | for modeling absorption effects. However, our limits do 
not depend much on the adopted background light modelQ since th ey a re dominated by the 



lOgl p. When calculating 



Galactic signal in most cases (for details of the calculation see Ref. 
the ICS component, however, we conservatively only include the extragalactic part, coming from 
electrons/positrons from dark matter decay that scatters on the CMB. The calculation of the 
ICS emission inside our Galactic diffusion zone is plagued with uncertainties and a detailed 



study is beyond the scope of this paper (see e.g. Ref. [109] for a thorough discussion). When 
quoting limits, we will require that in none of the energy bins considered in Ref. [lOll ] the dark 
matter signal integrated over these bins exceeds the m easur ed flux by more than 2a. Such limits 

or subt racting known astrophysical 

contributions to the extragalactic gamma-ray background 



llll.lll2l|. 



4 Results 

In none of the galaxy clusters a gamma-ray emission was found at the 3cr level, neither when 
searching for decay nor for annihilation signals with or without substructure contributions El We 
derived 95% C.L. limits on the dark matter lifetime and annihilation cross-section, from individual 
clusters as well as in a combined likelihood analysis; our results are shown in Figs. HI [U [6] and [71 
assuming 100% branching ratios into 66, fi~^fi~, W~^W~ or r+r" final states. 

Main Results. For decay or annihilation into 66, as relevant for MSSM neutralino DM, the 
strongest individual limits come from the Fornax cluster in most cases, as shown in Fig. [H (the 
impact of dark matter substructures on the annihilation limits is discussed below). Depending 
on the dark matter mass, lifetimes up to 4 • 10^^ s and annihilation cross-sections down to 5 • 
10~^^ cm^ s~^ can be constrained. Further strong limits come from AWM7, S636 and NGC4636. 



Our limits on the dark matter lifetime are somewhat weaker than previous results [27|; the 
difference can be mainly attributed to the fact that we modeled the cluster emission as an extended 
signal rather than as a point-like source, as we will discuss below. 

The limits obtained from our combined likelihood analysis are shown as dashed black lines in 
Fig. [3J They are often slightly weaker than the strongest individual limits. This is due to a weak 
preference for a non-zero signal in some of the clusters (like in A1367). In any case the combined 



^■^The adopted background light model appears to be in conflict with recent Fermi LAT observations, see 
Ref. [Tost . 

^•^The best signal candidate comes from A1367 (annihilation into r+r^, jtidm ~ 10 GeV) with a trial-corrected 
significance of 2.7a. 
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Figure 4: Left panel: lower limits on dark matter lifetime for decay into bb final states, as function of 
the dark matter mass. Solid lines show individual cluster limits, the dashed line the limits from the 
combined likelihood analysis. The dotted line shows for comparison the limit that can be derived from 
the EGBG. Riglit panel: like left panel, but upper limits on annihilation cross-section. Contributions 
to the signal from dark matter substructures are neglected. Note that the combined limits are 
sometimes weaker than the strongest individual limits; this is due to a weak preference of a non-zero 
signal for some of the clusters, see discussion in text. 



likelihood limit is more robust with respect to uncertainties of the cluster masses, the background 
modeling and statistical fluctuations in the data than the individual limits (see below discussion 
and Fig. E]). 

For comparison, the dotted line in the left panel of Fig. H] shows the lifetime limits that we 
obtain from conservatively requiring that the isotropic dark matter signal does not overshoot the 



EGBG as determined by Fermi LAT 10l[ |. see Section [3l The EGBG limit clearly dominates the 
cluster limits for large dark matter masses, whereas the cluster limits are competitive for masses 
below a few hundred GeV. In any case, since the systematics related to background subtractions 
are different for EGBG and cluster limits, the limits should be considered as being complementary. 

Limits on decay or annihilation into /x"//"*" final states are shown in Fig. [5l This channel is 
relevant for leptophilic models 57-63|, that aim to explain the PAMELA/Fermi anomalies. 



the corresponding best fit regions being shown in green (PAMELA only) and blue (PAMELA + 
Fermi + H.E.S.S.) In the presented dark matter mass range, the dark matter signal is 

dominated by ICS radiation of the produced electrons and positrons on the CMB; the prompt 
final-state-radiation can be neglected. Our dark matter lifetime limits reach up to 3 • 10^^ s for 
individual clusters as well as in the combined likelihood analysis, with the strongest limit coming 
from Fornax. In the case of dark matter annihilation limits down to 6 x 10"^'^ cm'^ s~^ are obtained. 
The parameter space favored by PAMELA /Fermi is constrained but not excluded in case of dark 
matter decay, and remains practically unconstrained in case of dark matter annihilation. 

The dotted line in the left panel of Fig. [5]shows again the limit obtained from the EGBG. In the 
case of decay into the cluster lifetime limits actually dominate over our conservative EGBG 

limit at all considered dark matter masses. This is due to the fact that we neglected the Galactic 



^■^Recently released Fermi LAT results [ll3l | indicate that the positron fraction continues to rise up to energies 
of 200 GeV, which will presumably shift the prefered DM mass range to somewhat higher values. 
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Lower limits on DM decay rate; 95% C.L. 



Upper limits on DM annihilation rate; 95% C.L. 




Figure 5: Like Fig. [4l but for jJ.^ fJ.'^ Rnal states. The blue region indicates the parameter region 
where a good fit to the PAMELA/Fermi LAT and H.E.S.S. electron/positron data is achieved \6aj . 
the green region corresponds to a fit to the PAMELA data only. 



ICS emission when calculating the EGBG limit. Such a calculation would require a treatment of 
cosmic-ray propagation in our Galax y, which has its own specific uncertainties and is beyond the 
scope of this paper (see e.g. Ref. |inj v On the other hand, as discussed above, propagation effects 
in galaxy clusters can be neglected at angular scales relevant for LAT observations, making the 
cluster limits on /i"*"^" final states practically independent of propagation model uncertainties. 




Figure 6: Left panel: lower limits on dark matter decay rate from the combined likelihood analysis, 
for different Bnal states, as function of the dark matter mass. We show limits at 95% C.L. (solid) 
and 99.7% C.L. (dotted) for comparison. Right panel: corresponding upper limits on annihilation 
rate. Contributions to the signal from dark matter substructures are neglected in this plot. 

In Fig. [6l we finally present a summary of our combined likelihood limits on dark matter 
decay and annihilation into different final states, bb, W~^W~ and t~^t~ . For comparison, we plot 
the 95% C.L. as well as the 99.7% C.L. limits. Limits on W^W~ are up to mass-independent 
rescaling similar to the limits on 66; the limits on t'^t~ clearly indicate that at dark matter 
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masses of ~ 1 TeV the ICS part of the dark matter signal starts to dominate inside the considered 
gamma-ray energy range. 




Figure 7: Left panel: individual upper limits on annihilation rate into bb, including a signal boost 
from dark matter substructures as discussed in Section \2. 31 (solid lines). For comparison, the dotted 
lines show the limits on the unboosted signal, cp. Fig. [4] Right panel: the same, but for annihilation 
into T^T~ . 



Dark Matter Substructures. In the above hmits we neglected contributions from dark matter 
substructures to the dark matter annihilation signal, firstly to obtain very conservative limits 
and secondly for the sake of comparison with previous work. However, dark matter substructures 
are a prediction of cold dark matter scenarios and expected to boost the annihilation signal 
considerably (29I . 30, 93, 94, [96|. To study their possible impact on our limits, we follow the 



prescription presented in Ref. [30|, which builds on results from the Aquarius project j92l . |93| 
and leaves the free streaming mass scale as the only free parameter (see Sec. 12.31 for a 

discussion). The resulting signal profiles are plotted in the right panel of Fig. [Uby the dotted 
lines, where we adopted a free streaming mass of Miim = W-^Mq. As can be seen from this plot, 
the boosted signal profiles extend to much larger radii than the profiles from the smooth dark 
matter halo alone. The corresponding boost factors of the overall signal are of the order of 10' 



,3 



consistent with what is found in Refs. [3d. I96l|. We note, however, that in the recent literature 



also smaller boost factors for galaxy clusters were discussed @, [23] (see also discussion above) , 
and our adopted boosted fiuxes should be considered as being optimistic but not unrealistic. 

No evidence for an extended annihilation signal due to dark matter substructures was found. 
In Fig. [7] we show the corresponding 95% C.L. limits on the annihilation cross-section compared 
to the limits obtained without dark matter substructures taken into account (for simplicity we 
neglect uncertainties in the overall cluster mass in case of the boosted signal). As expected, 
we find that the limits are improved by a factor of a few hundred; in the case of dark matter 
annihilating into bb with a thermal cross-section of 3 x 10~^^ cm^ s~^, dark matter masses below 
~ 150 GeV can be excluded. For different values of the limits would approximately scale 
like oc M^^^. 
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Lower limits on DM decay; Comparison 



Upper limits on DM annihilation; Comparison 




Figure 8: Like Fig. [4l but comparison of some of our results (solid lines, cp. Fig. \^ with the limits 
we would obtain when neglecting cluster mass uncertainties (dashed lines), and when furthermore 
approximating the cluster emission by a point-like source (dotted lines). 



Discussion. Comparing our results to previously presented limits 2a, 27|, we find that despite 
the increased statistics (at least by a factor of three) our limits on dark matter decay into bb 



remain somewhat weaker than the dark matter lifetime limits presented in Ref. [271 ]. This is due 



to the point-source approximation of the dark matter signal which was underlying the analysis 



of Ref. [27l |. as well as the inclusion of cluster mass uncertainties in the present work. In the case 
of dark matter annihilation our limits improve the results from Ref. [2^ by up to a factor of two. 

In Fig. [8] we compare our 95% C.L. limits on decay or annihilation into bb (solid lines) with the 
limits we would obtain when dropping the cluster mass uncertainties (dashed lines), and when 
furthermore approximating the dark matter signal as a point-like source (dotted lines). In the 
latter case, the extended signals are integrated over a re gion with a radius of 1° and modeled 



as a point source at the cluster center, following Ref. [26l . |27| |. As can be seen from Fig. [8l the 
non-inclusion of uncertainties in the cluster masses could strengthen the limits by up to ~ 50% in 
some cases, similar for both dark matter decay and annihilation. A point-source approximation to 
the signal from decaying dark matter would lead to a further strengthening of the corresponding 
limits by a factor of a few in some cases, as shown in the left panel of Fig. [8] (a similar impact is 
expected in case of annihilation signals boosted by dark matter substructures). However, in case 
of an unboosted annihilation signal (right panel), the point-source approximation to the signal 
changes the limits only on the 10%-30% level and hence appears to be justified, in agreement 
with what was found in Ref. \2i 



The robustness of a limit with respect to the underlying statistics and uncertainties can be 
inferred from its dependence on the adopted C.L. In Fig. [9] we show for the case of dark matter 
decay into bb our individual and combined limits at 95%, 99.7% and 99.994% C.L. in comparison. 
As can be seen from these plots, the combined limit depends considerably less on the adopted 
C.L. than the individual limits. This is due to the fact that uncertainties related to the cluster 
mass, as well as statistical fluctuations in the target region, tend to average out in the combined 
likelihood analysis. Furthermore, note that systematic effects in the modeling of the astrophysical 
backgrounds become in general less important at higher C.L. 
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Lower limits on DM decay rate at different C.L. 
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Figure 9: Like left panel of Fig.\^ but at different C.L. Note that the combined limit is less dependent 
on the C.L. than the individual limits, since uncertainties in the cluster masses tend to average out 
if limits are combined. 



In this paper we adopted galaxy cluster masses that are based on the X-ray observations from 



Refs. 69|, |70(], the HIFLUGCS catalog. The main advantage of this catal og is its c ompleteness. 



Many other X-ray measurements of cluster masses exist (see e.g. Refs. 114 - lid ]), which are 



however often concentrating on only one single cluster at a time, and hence difficult to use 
in a combined analysis. Important assumptions that enter the mass determination via X-ray 
observations is that the intra-cluster gas is in hydrostatic equilibrium and spherically distributed; 
further assumptions concern the temperature gradient of the gas, which is often approximated 
to be zero, as well as the radial distribution of the gas density. The systematical errors made by 
these approximations are very difficult to estimate, and it is likely that a neglect leads to a too 
optimistic determination of the cluster mass uncertainties. 

Fortunately, X-ray measurements are not the only way to determine the mass of a galaxy 
cluster; other methods include studies of the veloci ty dispersion of cluster member galaxies and 
weak gravitational lensing (see e.g. Refs. 

A comparison of the masses derived from 



different methods can be used as a rough estimate for the overall systematic errors in the cluster 



mass measurements. In Refs. j26l | and [271], the dark matter signal fluxes that follow from different 
cluster mass measurements (namely for the M49, Fornax and Coma clusters) were compared, and 
it was concluded that the overall uncertainties are roughly of the order of a factor of ~ 2 and 
not systematically biased with respect to the HIFLUGCS X-ray values. Looking at Tab. [U one 
can see that this is somewhat larger than what follows from the mass uncertainties given in the 
HIFLUGCS catalog; the latter are derived from uncertainties in X-ray profile fit, the temperature 
measurements and the temperature gradient. Hence, we expect that our limits on the dark 
matter annihilation cross section and lifetimes are not the most conservative ones that one could 
obtain for individual clusters when exploiting all systematic uncertainties. The advantage of our 
combined analysis is that, in absence of a systematic bias of the HIFLUGCS catalog, these kind 
of uncertainties are expected to partially average out. 

In the above calculation of the dark matter induced ICS emission the possible impact of 
intra-cluster magnetic fields was neglected; this is only justified if the magnetic fields remain well 
below the critical value of i?cMB = 3.2/iG. However, since Faraday rotation based measurements 
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in galaxy clusters find magnetic fields of a few fiG (see e.g. 0, 121 . 1221 ]). the validity of this 
approximation is not guaranteed. For most of our targets in Tab. [T]the magnetic fields are not 
precisely know, which makes it difficult to systematically include their effec ts. However, in case 
of the very massive Coma cluster the magnetic field was studied in Ref. 1211 ] , and a model for the 
magnetic field profile was presented. Adopting this model we can calculate how our limits change 
when synchrotron emission of the electrons and positrons produced in the dark matter decay or 
annihilation is included. We find that in case of the Coma cluster the limits on dark matter 
annihilation into as shown in Fig. [5] are weakened by a factor of around two, whereas the 

effect on decaying dark matter limits is negligible (the same is true for the extended annihilation 
signals coming from dark matter substructures). In the adopted magnetic field model, the field 
exceeds the critical value i?cMB only close to the cluster center at angles 6 < 0.17°, making 
the impact on point-source like signals large and on extended signals small. We conclude that 
intra-cluster magnetic fields are unlikely to affect our decaying dark matter limits or the limits 
on subhalo-boosted annihilation signals. 



5 Consequences for gravitino dark matter 



The gravitino is the spin-3/2 supersymmetric partn er of the graviton. If the lightest superparticle 
(LSP), it provides a natural dark matter candidate 123| . the mass of which can vary from 7113/2 ~ 
eV to ~ TeV depending on the details of the supersymmetry breaking mechanism. Gravitinos 
are produced in the early universe through 2-to-2 thermal scatterings with an abundance which 
is proportional to the reheating temperature Tr after infiation 



C 



100 GeV^ 



m 



3/2 



rrin 



[ Tr 



ITeV/ VlQiOGeV 



(14) 



where m^i2 and rrig are the gravitino and gluino masses respectively, and C ~ 0.5 to leading 
order in the gauge couplings 124 - 126 ]^ In addition, gravitinos may also be produced through 
the gravitational decay o f the NLSP. However, for r^NLSP^^ <C 1 or mNLSP ^ ^^3/2 the latter 
contribution is n egligible [1271 ]. Moreover, infiaton decay may also contribute to the production 



mechanism 12^ . In what follows, thermal leptogenesis is assumed to be responsible for the 
generation of the observed baryon asymmetry. In such a case, high reheating temperatures are 
required, and th e domina nt gravitino production mechanism is the thermal one. Particularly, for 
Tr ~ lO^o GeV lii^,^ a gravitino abundance of the order of the observed dark matter relic 
density Q-dm = 0.11 [tJis achieved for typical supersymmetric parameters, i.e. m3/2 ~ 100 GeV 
and rrig ~ 1 TeV. However, as it is well known, such high values of the gravitino mass lead to 
slow NLSP decays and can dram atically a ffect the successful predictions of the standard big bang 
nucleosynthesis (BBN) scenario [13ll4l37l ]. 

Among the different scenarios proposed to reconcile thermal leptogenesis, gravitino dark mat- 
ter and BBN, a mild violation of i?-parity inducing a rapid decay of the NLSP before the onset of 
the BBN is of interest [31|. In such a case, the gravitino is not stable anymore, but still provides 
a viable dark matter candidate due to the double suppression of its decay, by the Planck scale 
as well as by the small i?-parity breaking parameter. Interestingly, this opens up the way to 
look for traces of gravitino decays in cosmic-ray fiuxes, such as anti-matter [36|, |37[ and neutrino 
[3^. Additionally to the intense gamma-ray line arising from the 'ip?,/2 ~^ 7^ two-body decay 



^Note that C has 0(1) uncertainty due to unknown higher order contributions and nonperturbative effects [12 
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3ll435l. l37l ]. the produced gamma-ray flux typically features a continuous component generated 
by the fragmentation of the Higgs and gauge bosons. 

In what follows, we apply the above analysis to the decaying gravitino scenario. Contrarily to 
gamma-ray lines, galaxy clusters offer more sensitivity to large gravitino masses, thus rendering 



the present analysis supplementary to our previous gamma-ray lines study [47[. Following the 
structure of the latter, we first summarize the bilinear i?-parity violation supersymmetric frame- 
work considered here. We then present limits on the size of i2-parity violation and finally discuss 
the prospect for seeing long-lived neutralino and stau NLSPs at the LHC. 

5.1 i?-parity breaking model 

The supersymmetric standard model with explicit bilinear i?-parity violation is specified by the 
superpotential 

W = Wmssm + ^J'^Huli , (15) 
as well as by the soft supersymmetry breaking potential 

C = + B.Huk + mjJlH^ + h.c. , (16) 

where VFmssm and C^^^^^^ are the i?-parity conserving MSSM superpotential and scalar La- 



grangian, i/^/^ are the up/down- type Higgs doublets, k the lepton doublets, and fii, Bi and 
are the -R-parity violating couplings. Trading the mass mixing parameters for i?-parity breaking 



id 



Yukawa couplings as proposed in Ref. 1381 ] . the gravitino decay is function of a sin gle d imen- 



sionless parameter C,, which also enters the decay of the NLSPs of interest (see Ref. 138l | for a 



definition of Q in terms of the bilinear i?-parity violating couplings /Xj, Bi and mnf^. 

Two typical sets of boundary conditions for the supersymmetry breaking parameters of the 
MSSM at the grand unification (GUT) scale are investigated in the following, resulting in two 
different types of NLSPs. First we consider equal scalar and gaugino masses 

(A) mo = mi/2) «o = 0, tan/3 = 10, (17) 

for which the bino-like neutralino is the NLSP. In the second one, which corresponds to no-scale 
models or gaugino mediation, 

(B) mo = 0, mi/2/0, ao = 0, tan/3 = 10, (18) 

the lightest stau fi is the NLSP. In both cases, tan/3 = 10 has been chosen as a representative 
value, and the trilinear scalar coupling oq has been set to zero for simplicity. For both sets of 
boundary conditions, the universal gaugino mass mi 12 remains as the only independent variable, 
and the gaugino masses Mi^2,3 satisfy the following relations at the electroweak scale 

M3 M2 

Electr oweak precision tests (EWPT) yield important lower bounds on the superparticle mass 



spectrum |126l |. For a neutralino NLSP, the most stringent constraint comes from the Higgs 



potential. The universal gaugino mass mi/2 is required to be high enough in order for the 



Higgs mass to fulfills the LEP lower bound rrih > 114.4 GeV jl39l |. This implies the lower limit 



m^o > 130 GeV cj However, allowing negative oq or scalar masses much larger than mi/2 at the 



16 



Note that m^o ~ Mi with good accuracy [13i 
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GUT scale would weaken this limit, and we will take m^o > 100 GeV as a lower bound for the 
neutralino mass subsequently. In the stau NLSP case, the lower bound comes f rom the absence 



of pair production of heavy charged particles at LEP and reads rrif-^ > 100 GeV [139l |. Rewriting 
Eq. (mi) 



where ^ = rriNLSP /"T-g is implicitly fixed by the supersymmetry breaking boundary conditions 12d | , 
we get absolute upper bounds on the NLSP masses requiring the gravitino to be the LSP. In the 
case of the neutralino NLSP, Eq. (pO]) implies m^o < 690 GeV for ^ = 1/5.9, and is essentially 
independent of tuq and tan /3. For the stau NLSP, tan /3 = 10 yields ^ = 1/6.2, which consequently 
leads to the more strin gent bounds < 615 GeV. Note that there is a strong dependence on 



tan/3 in that case [126l |. and that ^ decreases with increasing tan/?. 

For a typical effective neutrino mass fhi = 1Q~^ eV, successful thermal leptogenesis requires 



a minimal reheating temperature of ~ lO^GeV H- Using Eq. (1140 together with a lower 



bound on the gluino mass nig > 815 GeV (l40j . this implies a lower bound for the gravitino mass 



m 



3/2 



> 30 GeV. 



5.2 Limits from galaxy clusters 

The gamma-ray spectrum produced through gravitino decays features two types of contributions: 
First, the ~^ 1^ two body decay produces a gamma-ray line, a channel which is dominant 
below the W thresholc0. Additionally, both the fragmentation of the Higgs and gauge bosons as 
in Z^v and h^v final states, as well as the final state radiation of the charged leptons produced in 
the W^r^ final states, generate a continuum spectrum. The relative strength of these two is fixed 
by the corresponding branching ratios, that we present in Fig. [10] together with representative 
spectra following Refs. 37|, |39||^ While the search for gamma-ray lines of galactic origin thr ough 



deviations from a power law background is efficient for contained continuum contributions [47|, 
i.e. for m3/2 ^ 200 — 300 GeV, constraints from galaxy clusters observations and the EGBG 
dominates the gamma-ray line ones for 777-3/2 ^ 250 GeV, as illustrated in Fig. [TTJ This agreeably 
makes gamma-ray line searches, galaxy cluster observations and EGBG studies complementary. 
As a result the gravitino lifetime is constrained to be at least r3/2 ^ 0(10^^ s) in all the gravitino 
mass range considered. Subsequently, we will concentrate on the limits derived by galaxy cluster 
observations, and leave a detailed study of implications from the EGBG to future work. 

5.2.1 i?-parity breaking parameter 



The gravitino inverse decay rate into photon/neutrino pairs is given by 321. Il38l | 



, 32^/2 GpMg MfM| / / m|\\ 

where a is the electromagnetic fine structure constant, Mp = 2.4 x 10^^ GeV the reduced Planck 
mass, and Gf = 1.16 x 10"^ GeV~^ is the Fermi constant. Using the strongest limits on the total 



^^Note that three-body decays with inter medi ate massive gauge bosons are expected to contribute by more than 
10% below the kinematic threshold [i^. lii. [l4lt . 

^*Note that the branching ratio into lines is in principle model-dependent. 
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Figure 10: Left: Two-body decay branching ratios of the gravitino. Right: Gamma-ray spectra for 
"^3/2 = 100, 200, 500 and 1000 GeV. We adopt here the same set of parameters as in Refs. 0, [g^j . 

gravitino lifetime illustrated in Fig. [TTjtogether with the branching ratios presented in Fig. [TUl this 
expression can be used to derive conservative upper-limits on the i?-parity breaking parameter 
^. To do so, one has to consider for a given gravitino mass the maximally allowed bino mass 
which results from the combination of Eqs. (jl4p and (jl9p when considering the lowest reheating 
temperature allowed in the thermal leptogenesis scenario, i.e. Tr ~ 10^ GeV. The results are 
presented in Fig. [T2j Note that at high gravitino masses, the production of anti-protons in Higgs 
and gauge bosons fragmentation could further constrain the ^ parameter. 

5.2.2 Stau NLSP decay length 

In the case of the fi-NLSP, the total decay width of the lightest mass eigenstate is a mixture of 
left and right handed partial decays 

Tr, (e) = sin^ 9^ Tf^ (e) + cos^ Of Tf^ (e) . (22) 

where the dim ensionless parameter e is directly related to the i?-parity violating Yukawa couplings 



(see Ref. 138l | for details). Since the latter are typically proportional to the ordinary Yukawa 



couplings, decays into second and third families dominate. For definiteness, we will below assume 



a flavor structure as described in Ref. 138l |. into which the chiral state decays are dominated by 
the following channels 

TR TLf, ULiy , (23a) 
TL ^ tnbL . (23b) 

Assuming ~ e0 and using the cluster upper limits on ( from Fig. \T2\ we can derive lower 
bounds on the stau decay length. Our results are shown in Fig. [131 The parameter space is already 



% values much smaller than e can be achieved through a proper choice ol the parameters /li, Bi and rrq^. 



22 




'3/2 



[GeV] 



Figure 11: Lower limits on tlie gravitino lifetime. The dot-dashed line shows the gamma-ray line 
limits, the dashed line the limits resulting from the combined cluster analysis, and the dotted line 
the EGBG limits. 



constrained by EWPT and overproduction bounds, and the lower limits on the neutralino decay 
length vary between 100 m and 10 km. It is interesting that if such particles wer e to be produced 
at the LHC, a sizable amount of their decays could take place in the detector |38l . ll42l |. We obtain 
a lowest possible decay length cTf-^ ~ 200 m for 



30 GeV and nif-^ 



130 GeV. 



5.2.3 Neutralino NLSP decay length 



A neutralino NLSP heavier than 100 GeV dominantly decays into W^t^ and Z^u 143l . ll44l ]. The 
corresponding decay width is directly proportional to the i2-parity breaking parameter <^ squared, 
which also enters t he gr avitino decay width Eq. (i2T]l . As a consequence, the two quantities can 
be related through |l38l | 



cl {M2-Mif 



m' 



3/2 



r 



2^/2 



G-FM^mlo 2f{m^o,mw) + f{'m^o,mz) 



(24) 



where the phase space factor / is defined by 

f {mi, 1712) = 



ml 



1 + 24 

mf 



(25) 



Using the gaugino mass relation Eq. (|19p . lower bounds on the neutralino decay length cr^o 
can be derived from the partial gravitino decay width. Our results are summarized in Fig. [ril 
considering the cluster lifetime limits. For the parameter space allowed by EWPT and overpro- 
duction bounds, we obtain mini mal decay lengths O(100m — 100km), which are in the range 



of detectability of the LHC [38|, Il42l |. 

7713/2 ^ 30 GeV at m^o ~ 140 GeV. 



Decay lengths as small as cr^o 



60 m are allowed for 
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Figure 12; Upper bounds on the R-parity violation parameter derived from the Fermi LAT hmits 
in Fig. 11 Jl For thermal leptogenesis, overproduction (O.P.) of gravitinos already excludes the left 
green region. The dotted line illustrates the limits from EGBG. 

6 Conclusions 

Galaxy clusters are the most massive collapsed objects in the Universe, and very interesting tar- 
gets for the indirect search for dark matter. Using three years of Fermi LAT data, we studied 
the gamma-ray emission from eight of the most promising galaxy clusters and searched for signa- 
tures for dark matter decay or annihilation. We analyzed the clusters individually as well as in 
a combined likelihood approach. We took into account uncertainties in the cluster masses as de- 
termined by ROSAT PSPC X-ray observations and modeled the dark matter signals as extended 
sources. Our main results are: 

• In none of the eight galaxy clusters listed in Tab. [J a significant gamma-ray emission that 
could be attributed to dark matter decay or annihilation was found. We derived limits 
at the 95% C.L. on the dark matter lifetime and on the annihilation cross-section, from 
each cluster individually as well as in a combined likelihood approach. Our results are 
shown in the Figs. HI [5l [6] and [7] for bb, W'^W~ and T~'~r~ final states. In most 
cases the combined limits are at the level of the strongest individual limits; in any case 
the combined limits are more robust with respect to uncertainties of the cluster masses, 
background modeling and statistical fluctuations in the cluster target regions (Fig. [9]). 

• The limits on the dark matter lifetime turn out to be somewhat weaker than what previous 
results indicated 271], reaching up to lifetimes of tdm — 4 x 10^^ s in case of decay into bb 



(Fig. H]). The difference can be partly attributed to the fact that we modeled the cluster 
emission as an extended signal rather than as a point-like source, and partly to the inclusion 
of cluster mass uncertainties. We find that a point-source approximation to the signal from 
dark matter decay could in some cases strengthen the limits wrongly by a factor of a few 
(Fig. [8]). In particular, the decaying dark matter interpretation of the excess in terms 
of dark matter decaying into /i'''/^" remains only partially constraint by our galaxy cluster 
limits (Fig. [5|). 

As long as prompt radiation dominates the overall gamma-ray signal, our cluster limits 
on the dark matter lifetime are sometimes weaker than the corresponding limits that can 
be conservatively inferred from measurements of the extragalactic gamma-ray background 
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Figure 13: Contour plot of lower bounds on the stau NLSP decay length coming from cluster and 
gamma-ray line constraints on the gravitino lifetime, as function of the stau and gravitino masses, mf^ 
and m3/2 respectively. The lower gray region is excluded by electroweak precision tests (EWPT). For 
thermal leptogenesis, overproduction (O.P.) of gravitinos excludes at minimum the left green region, 
a limit which strengthens to the black-dashed line when assuming the universal boundary conditions 
(B), cf. Eq. ((TSJ. The lower-right exclusion comes from the gravitino LSP requirement. 



(EGBG, see Figs. H]and[5|). However, at lower gamma-ray energies the cluster and EGBG 
limits are comparable and should be considered as being complementary, since the system- 
atics of the background subtraction are in general very different. 

• In case of an annihilation signal coming from the smooth component of the cluster dark 
matter halo alone, while neglecting the contributions from dark matter substructures, pre- 



viously presented limits (26|] are improved by up to a factor of two (Figs. HI [5] and [6]). In case 
of annihilation into bb final states, we obtain limits down to 5 x 10~^^ cm^ s"'^ for a 10 GeV 
WIMP, and we confirm that the point-source approximation to the dark matter emission 
is valid when the effect of dark matter substructures are neglected. However, as we discuss 
for the case of the Coma cluster, for such a point-like emission intra-cluster magnetic fields 
could reduce the unboosted annihilation signal by an order one factor when the signal is 
dominated by photons from ICS radiation. 

• When contributions to the annihilation signal from dark matter substructures are taken 
into account, the dark matter signal must be modeled as an extended source {cp. Fig . [T]). 
Adopting the optimistic scenario for signal boosts due to substructures from Ref. [SO], we 
find that the limits could strengthen by a factor of a few hundred if substructures with 
masses down to Afum = 10~^Mq are included. In this case, the limits on bb would start to 
reach the thermal cross-section 3 x 10~^^ cm^ s~^ and exclude dark matter masses below 
150 GeV (Fig. ED. 

As a direct application of our results, we derived limits on the decaying gravitino dark matter 
scenario, both from galaxy clusters observations and from the EGBG. We find that the cluster 
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Figure 14: Like Fig. \T3[ but for a neutralino NLSP. 

limits on the gravitino lifetime is ^ x 10^^ s for gravitino masses up to 1 TeV. The cluster 
constraint becomes stronger than the gamma-ray line limit at gravitino masses of 7713/2 ^ 250 GeV. 
As a result, we found limits on the ii-parity breaking parameter of the order of 0(10"^*') for 
gravitino masses 250GeV< 7713/2 ^ 1 TeV. These constraints were used to set lower limits on 
NLSP decay lengths corresponding to two different types of universal boundary conditions at the 
grand unification scale of supergravity models. Interestingly, all the implied decay lengths should 
be accessible at the LHC. 
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